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The self-force describes the effect of a particle's own gravitational field on its motion. While the 
motion is geodesic in the test-mass limit, it is accelerated to first order in the particle's mass. In this 
contribution I review the foundations of the self-force, and show how the motion of a small black 
hole can be determined by matched asymptotic expansions of a perturbed metric. I next consider 
the case of a point mass, and show that while the retarded field is singular on the world line, it can 
be unambiguously decomposed into a singular piece that exerts no force, and a smooth remainder 
that is responsible for the acceleration. I also describe the recent efforts, by a number of workers, 
to compute the self-force in the case of a small body moving in the field of a much more massive 
black hole. 

The motivation for this work is provided in part by the Laser Interferometer Space Antenna, 
which will be sensitive to low- frequency gravitational waves. Among the sources for this detector is 
the motion of small compact objects around massive (galactic) black holes. To calculate the waves 
emitted by such systems requires a detailed understanding of the motion, beyond the test-mass 
approximation. 

This article is based on a plenary lecture presented at the 17th International Conference on 
General Relativity and Gravitation, which took place in July, 2004 in Dublin, Ireland. 



INTRODUCTION: THE CAPRA SCIENTIFIC 
MANDATE 



This contribution describes how a body of mass m, 
supposed to be "small" , moves in spacetime. In the test- 
mass approximation it is known that the body moves 
on a geodesic of a background spacetime whose metric 
g a fj does not depend on m. As m increases (but kept 
small) the test-mass description is no longer adequate. 
One might then say that the motion is still geodesic, but 
in a spacetime whose metric is perturbed with respect to 
the background metric. Alternatively, one might say that 
the motion is accelerated in the background spacetime. 
This is the point of view I shall take in this contribution: 
The body's motion will be described in the original (un- 
perturbed) spacetime, and the gravitational influence of 
the body will be incorporated in its acceleration. The 
body will be said to move under the influence of its grav- 
itational self-force. 

The work described here was previously reviewed in 
a very long article published in Living Reviews in Rel- 
ativity [l| . An abridged version of this article appeared 
in Classical and Quantum Gravity 2]. This contribution 
focuses on the highlights and presents the "big picture" . 
The reader is warned that my presentation will be sketchy 
and incomplete, and is referred to the review articles for 
additional details. The scientific objectives behind the 
work described here have been pursued by a number of 
people; I like to refer to these objectives as the "Capra 
scientific mandate" and to these people as the "Capra 
posse" . 

The Capra scientific mandate includes three main ob- 
jectives. The first is to formulate the equations of mo- 
tion of a small body in a specified background spacetime, 
beyond the test-mass approximation. This first step was 
solved in 1997, first by Mino, Sasaki, and Tanaka j3j, and 



then by Quinn and Wald . The equations of motion are 
now known as the MiSaTaQuWa equations; I will sketch 
their derivation in this contribution. The second objec- 
tive is to concretely describe the motion of the small body 
in situations of astrophysical interest, including generic 
orbits of a Kerr black hole. There has been much re- 
cent progress on this front, and I will describe some of 
the issues involved in this contribution. The third ob- 
jective is to properly incorporate the equations of motion 
into a wave-generation formalism. This final objective 
is the most important, as the ultimate goal of this en- 
terprise is to make detailed predictions toward eventual 
gravitational-wave measurements. This is the holy grail 
of the Capra program, and it has so far proved elusive. I 
will present a tentative outline of future work in the last 
section of this contribution. 

The work reviewed in this contribution was shaped 
by a series of annual meetings named after the late 
movie director Frank Capra. The first of these meet- 
ings took place in 1998 and was held at Capra's ranch in 
Southern California; the ranch now belongs to Caltech, 
Capra's alma mater. Subsequent meetings were held in 
Dublin, Pasadena, Potsdam, State College PA, Kyoto, 
and Brownsville. 

Members of the Capra posse include Paul Anderson, 
Warren Anderson, Leor Barack, Patrick Brady, Lior 
Burko, Manuella Campanelli, Steve Detweiler, Eanna 
Flanagan, Costas Glempedakis, Abraham Harte, Wataru 
Hikida, Bei Lok Hu, Scott Hughes, Sanjay Jhingan, 
Dong-Hoon Kim, Carlos Lousto, Eirini Messaritaki, Ya- 
sushi Mino, Hiroyuki Nakano, Amos Ori, Ted Quinn, 
Eran Rosenthal, Norichika Sago, Misao Sasaki, Takahiro 
Tanaka, Bob Wald, Bernard Whiting, and Alan Wise- 
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II. ASTROPHYSICAL CONTEXT 

The motivation for the Capra program comes largely 
from the fact that solar-mass compact bodies moving 
around massive black holes have been identified as one 
of the promising sources of gravitational waves for the 
space-based interferometric detector LISA (Laser Inter- 
ferometer Space Antenna). (The case for this identifica- 
tion is made in the contributions by Sterl Phinney and 
Sir Martin Rees.) These systems involve highly eccen- 
tric, nonequatorial, and relativistic orbits around rapidly 
rotating black holes. The waves produced by these orbits 
will be rich in information concerning the strongest grav- 
itational fields in the Universe, and this information will 
be extractable from the LISA data stream. The extrac- 
tion, however, will depend on sophisticated data-analysis 
strategies that will rely on a detailed and accurate mod- 
cling of the source. This modeling involves formulating 
the equations of motion for the small body in the field 
of the rotating black hole, in a small-mass-ratio approx- 
imation that goes beyond the test-mass description. And 
it involves a consistent incorporation of these equations 
of motion into a wave-generation formalism. In short, 
the extraction of this wealth of information relies on the 
successful completion of the Capra program. 

The finite-mass corrections to the orbiting body's mo- 
tion are important. Let m be the mass of the orbiting 
body, M the mass of the central black hole, and suppose 
that m/M <C 1. For concreteness, assume that the or- 
biting body is a 10 M & black hole and that the central 
black hole has a mass of 10 6 M©. Then m/M = 1CT 5 is 
the order of magnitude of the correction to the equations 
of motion relative to the test-mass description. Simul- 
taneously, M/m = 10 5 is the order of magnitude of the 
total number of wave cycles that will be received during 
a year's worth of LISA observation. This simplistic es- 
timate illustrates that while the corrections to the equa- 
tions of motion are small, in the course of a year they can 
accumulate and contribute a significant number of wave 
cycles. 

Corrections to the equations of motion must incorpo- 
rate both conservative and dissipative effects. Finite- 
mass corrections that are conservative in nature are fa- 
miliar from Newtonian and post-Newtonian theory, and 
they occur also in a strong-gravity situation. These can 
accumulate over time. Imagine, for example, an eccentric 
orbit that undergoes periastron advance. A finite-mass 
correction to this effect will cause a steady drift in the 
phasing of the orbit, and this will directly be reflected in 
the phasing of the gravitational waves; after 10 5 orbital 
cycles the 10 -5 correction will have grown into a sizable 
effect. 

Dissipative effects, on the other hand, do not occur 
in Newtonian gravity, but they are familiar in post- 
Newtonian theory; they are also present in a strong- 
gravity situation. Dissipation is associated with the ra- 
diative loss of energy and angular momentum by the or- 
bital system, and the resulting corrections to the motion 



also accumulate over time; this translates into a steady 
drift of the gravitational-wave signal in the LISA fre- 
quency band. These are finite-mass corrections, because 
the test-mass description makes no room for gravitational 
radiation and radiation reaction. 

From a more theoretical point of view, the appeal of 
this work comes largely from the fact that while the mo- 
tion of self-gravitating bodies has been studied exten- 
sively in the context of post-Newtonian theory (see, for 
example, Refs. 0,0); very little is known in the case of 
strong fields and fast motions. To the relativist work- 
ing in this area, this problem is irresistible: We have 
strong gravity, fast motion, and a smallness parameter 
in the form of m/M . We have a cool problem that can 
be solved by standard perturbative techniques. 

III. MOTION OF A BLACK HOLE 

Let me restate the problem in its most general form: 
A body of mass m moves in an arbitrary (but empty 
of matter) spacetime whose radius of curvature (in the 
body's neighbourhood) is 1Z; what is the description of its 
motion when m <C TZf This formulation of the problem 
is more general than the two-body version stated previ- 
ously. When the small body is a member of a binary 
system, and is at a distance b from another body of mass 
M, then K ~ y/b 3 /{M + m) and 

m m r> 

— ~ v 6 , 

11 M + m 

where v = a/ (M + m) Jb. For relativistic motion (v ~ 1) 
this is small whenever m/M <C 1. 

The clean separation of scales allows us to idealize the 
motion as following a world line in a spacetime whose 
metric g a p will be specified below. While the region of 
spacetime occupied by the body is truly a world tube of 
finite extension, on a scale 1Z — the only scale of rele- 
vance in the background spacetime — this extension is 
so small that little is lost by making this idealization. 
The body shall then follow a world line 7 that will be 
described by parametric relations z^(t), in which r is 
proper time measured in the metric g a p. We wish to de- 
termine this world line. It is understood that the world 
line loses all significance when the neighbourhood of the 
body is examined on the fine scale m; on this scale the 
finite extension of the world tube is fully revealed. 

It is desirable to choose the body to have an internal 
structure that is as simple as possible. One thus elimi- 
nates an important source of technical complications, but 
one still feels confident that the resulting equations of mo- 
tion will apply, to a good degree of accuracy, to a body of 
arbitrary structure. (This "effacement property" of the 
internal structure is well established in post-Newtonian 
theory |||.) This attitude often leads the researcher to 
assume that the body is a point particle. We shall re- 
frain from doing so at this stage, but we will come back 
to this description in Sec. IIVI We shall instead choose 
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the body to have the simplest structure compatible with 
the laws of general relativity: it shall be a nonrotating 
black hole. This was one of the starting assumptions of 
Mino, Sasaki, and Tanaka (see also Ref. 0). 

The motion of the black hole is determined by matched 
asymptotic expansions, a powerful technique that has 
known many useful applications in general relativity (see, 
for example, Refs. In this approach the metric of 

the black hole perturbed by the tidal gravitational field 
supplied by the external universe is matched to the metric 
of the external universe perturbed by the moving black 
hole. Demanding that this metric be a valid solution to 
the Einstein field equations determines the motion with- 
out additional input. 

The method of matched asymptotic expansions relies 
on the existence of (i) an internal zone in which grav- 
ity is dominated by the black hole, (ii) an external zone 
in which gravity is dominated by the conditions in the 
external universe, and (iii) an overlapping region known 
as the buffer zone, in which the black hole and the ex- 
ternal universe have comparable influence. Let r be a 
meaningful measure of spatial distance from the black 
hole. Then the internal zone is the region of spacetime 
in which r < = constant, where -C TZ; thus r/lZ 
is small throughout the internal zone and the tidal in- 
fluence of the external universe is small. The external 
zone, on the other hand, is the region of spacetime in 
which r > r e = constant, where r e 3> rn; thus m/r is 
small throughout the external zone and the gravitational 
effects of the black hole are small. Finally, the buffer 
zone is the region of spacetime in which r e < r < ri. 
This region exists because m C K, and r/TZ and m/r 
are simultaneously small in the buffer zone. 

An expansion of the metric in the internal zone has the 
schematic form 

g = <?(Schwarzschild) + 0(r 2 /K 2 ) + ■ ■ ■ , (1) 

where <?(Schwarzschild) is the metric of a nonrotating 
black hole in isolation, and 0(r 2 /7Z 2 ) is the tidal field 
supplied by the external universe. (Terms that scale as 
r/TZ have been removed by a coordinate transformation.) 
On the other hand, an expansion of the metric in the 
external zone has the schematic form 

g = g (background) + 0(m/r) + ■ ■ ■ , (2) 

where g (background) is the metric of the external uni- 
verse in the absence of the black hole, and 0(m/r) is the 
perturbation produced by the moving black hole. The 
key idea of matched asymptotic expansions is that the 
expansions of Eqs. (JIJ and 10) are both valid in the buffer 
zone, and that both forms of g represent the same met- 
ric (up to a coordinate transformation). Performing the 
matching returns the black hole's equations of motion. 



External zone 

To flesh out these ideas let us first examine the 
situation in the external zone. Let g a p stand for 
g (background), the metric of the external universe in the 
absence of the small black hole. (In the astrophysical 
context of Sec.^lthis would be the metric of the massive 
black hole, in isolation.) Let 7 be a fiducial world line 
in this background spacetime (later to be identified with 
the small hole's world line), and express the metric in 
normal coordinates centered on 7. It will have the form 
of g = r] + 0(r/lZ) + 0(r 2 /1Z 2 ) + - ■ ■, where r] is the metric 
of flat spacetime. The normal coordinates do not extend 
beyond r = TZ and are therefore restricted to the buffer 
zone. 

More specifically, let the normal coordinates be the re- 
tarded coordinates (u, x a — r£l a ), such that u is constant 
on each future light cone centered on 7 (and u is equal to 
proper time on the world line), r is an affine parameter 
on the cone's null generators, and il a = x a /r is constant 
on each generator. Then the time-time component of the 
background metric takes the form 

g uu = -(1 + 2ra a n a + r 2 £ ab n a n b ) + 0(r 3 /lZ 3 ), (3) 

where a a (u) is the acceleration vector of the world line, 
and £ a b(u) — C uau b — 0(1/1Z 2 ) are the electric compo- 
nents of the Weyl tensor evaluated on 7. (This three- 
tensor is symmetric and traccfrcc.) Our main goal is to 
eventually determine a a ; we can anticipate that a a = 
0(m/7Z 2 ) since the motion must be geodesic in the test- 
mass limit. 

Let h a p be the perturbation produced by the moving 
black hole, so that the full metric is g a p — g a p + h a p. A 
standard technique in perturbation theory is to introduce 
the trace-reversed potentials 

^ a = h a0 --g a0 (g' y5 h^s) (4) 
and to impose the Lorenz gauge condition 

= 0. (5) 

Here and below, all indices are manipulated with the 
background metric, and covariant differentiation is taken 
to be compatible with this metric. 

In the external zone the perturbation produced by the 
black hole cannot be distinguished from one produced by 
a point particle of mass m moving on 7, and ip a p must 
be a solution to the wave equation 

□V> Q/3 + 2R a y s ^ 5 = -167rT a/3 , (6) 

where T a @ is the particle's stress-energy tensor. This 
equation can be solved by means of a retarded Green's 
function, 

G a % s ,{x,x') = U a0 Y5l (x,x')S(<j)+V a ^ 5l (x,x')6(-<7), 

(7) 
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where a(x,x') is half the squared geodesic distance be- 
tween the points x and x', and U a ^, s , (x, x'), V a ^, 5 , (x, x') 
are smooth bitensors. The Green's function is decom- 
posed into a singular "light-cone part" that has support 
on a — only, and a smooth "tail part" that has sup- 
port on a < (so that x is in the chronological future of 
x'). (It should be noted that this representation of the 
Green's function is valid only when x is in the normal 
convex neighbourhood of x' . In the sequel I will simplify 
expressions by pretending that the decomposition holds 
globally. The reader is referred to LRR |]| for all the fine 
print.) 

The solution to Eq. (JBJ is 



4m . 



r v [x) = ^U a *, s ,(x,x>)uru* + (8) 

where x' = z^fu) is the point on the world line that 
is linked to a; by a null geodesic, u the value of the 
proper-time parameter at this retarded point, u 1 the 
four- velocity at the retarded point, and 

iCO*) = 4™ / V a P(x,z')u» u" dr' (9) 



with an extra factor of 1 /1Z, the inverse time scale asso- 
ciated with changes in the external universe. The pertur- 
bation must be well behaved on the hole's event horizon, 
and the asymptotic conditions for r ^> 2m are such that 
the perturbation must behave as — r 2 £ ab £l a £l b , which has 
a quadrupolar form. These observations make solving 
the perturbation equations a very straightforward task. 

In a set of internal light-cone coordinates (u, x a = fCt a ) 
the time-time component of the perturbed metric reads 



',m = -f- f 2 f 2 s ab n a n b + o(f 3 /n 3 ), 



(ii) 



where / = 1 — 2m jf and £ ab (u) is the tidal gravita- 
tional field supplied by the external universe. In the limit 
m/f — ► (keeping f/lZ fixed) Eq. becomes the met- 
ric of the external universe, expressed in coordinates for 
which a a = 0. In the limit f/lZ — > (keeping m/f fixed) 
Eq. (jllfl becomes the metric of a Schwarzschild black 
hole expressed in Eddington-Finkelstein coordinates. For 
small values of f/TZ and arbitrary values of m/f Eq. Hll|) 
describes a tidally distorted black hole. The metric of 
Eq. is valid everywhere in the internal zone, where 
r < r-j < ft. 



is the "tail" term. Here z' = z(t') stands for an arbi- 
trary position on the world line, and the potentials of 
Eq. (JSJ) depend on the particle's entire history prior to 
the retarded point r = u(x). 

Inverting Eq. Q and expressing the results in the re- 
tarded coordinates (w,rf2 a ) returns 



1 - r 2 £ ab n a n b + 0{r A /-R A ) 



2m 



*(2m£ ab n a tt b - 2a a n a 



+ h t u f u + h t u f a W)+0(mr 2 /ll 3 ) 



(10) 



IJ a/3 by trace reversal, and h a ^ 



The fields are obtained from 

tail 



Y7 h tail 



for the time-time component of the perturbed metric 

gaj3 — 9a(3 + h a p 
^tail 

n 

Equation l|10|) gives the metric of the external uni- 
verse perturbed by a black hole moving on a world line 
7 with an (as yet undetermined) acceleration a a (u). It is 
noteworthy that this metric appears to be singular when 
r — > 0. The metric, however, is valid only in the external 
zone r > 7' e 3> m, and the limit r — > is unattainable. 
The world line 7 is outside the metric's domain of valid- 
ity. 



Matching 

Equation (|10fl gives the spacetime metric in the exter- 
nal zone, and the requirement that r/TZ <C 1 implies that 
the metric is in fact restricted to the buffer zone. Equa- 
tion on the other hand, gives the spacetime metric 
in the internal zone, and its specialization to m/f < 1 
implies also a restriction to the buffer zone. Since both 
metrics describe the geometry of the same region of the 
same spacetime, they must be related by a coordinate 
transformation. 

The transformation from the external coordinates 
(u, x a ) to the internal coordinates (u, x a ) is given in part 
by 



u — 2m In r 



IL UU (ill 



r 



Ltail 

uuu 



1 

~2' 



■-.tail 



tail 

uua 



abu 



2hlf b )n a n b 



h t a f c n a n b n c 



+ (*2 

+ 0(mr 3 /K 3 ). 
Applying this transformation to Eq. 1|1(J[) gives 



(12) 



Internal zone 

To obtain the metric of a nonrotating black hole that is 
slightly distorted by a tidal gravitational field is a stan- 
dard application of black-hole perturbation theory. To 
leading order in the perturbation, which scales as 1/1Z 2 , 
it is sufficient to integrate a set of time-independent per- 
turbation equations, because each time derivative comes 



iuu = -i-f 2 £ ab Ci a Ci b + o(f 3 /n 3 ) 

2m 

+ — + Amf£ ab n a n b 
r 



-2r[a a 



1 



tail 



^tail 



+ 0{mf 2 /K 3 ). 



(13) 



Comparison with Eq. (|ll|l — linearized with respect to 
m/f — reveals that the acceleration vector must be given 
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by 

a a =^ a -hlf u . (14) 

As expected, the matching of the perturbed metrics in 
the buffer zone has returned the black hole's equations 
of motion. 



MiSaTaQuWa equations 

The tensorial form of Eq. I|14|) is 

= -~Gr + «"«") (?hX - h\f v )u x u», (15) 

and these are the MiSaTaQuWa equations of motion. 
Here, z^(t) gives the description of the black hole's mo- 
tion in the background spacetime, it M = dz^/dr is the 
normalized velocity vector, and Du^/dr is the covariant 
acceleration. All tensors, and all tensorial operations, re- 
fer to the background spacetime of the external universe; 
its metric is g a p, as it was introduced in Eq. ©. 

The tail field on the right-hand side of Eq. i|15|) is ob- 
tained by differentiating Eq. @ ; when expressed in terms 
of the original retarded Green's function of Eq. Q it is 
given by 

= 4m J r Jv A (<W"' -\g^G p pil , v )j{z,z!) 

x 1/1/ dT 1 . (16) 

It is understood that one takes the limit e — > + of this 
expression. Cutting the integration short excludes the 
(singular) light-cone part of the Green's function and iso- 
lates its tail part; the result is a smooth field on the world 
line. In Eq. (|16f) . z = z(t) stands for the current posi- 
tion on the world line (at which the self-force is being 
evaluated) and z' = z(t') stands for a prior position. 

IV. MOTION OF A POINT MASS 

The derivation of the MiSaTaQuWa equations 
sketched in the preceding section is conceptually sound 
and devoid of any questionable assumptions. It is, how- 
ever, technically involved and one wonders whether there 
may not be a shortcut to the final answer. Since the equa- 
tions of motion refer to a body without internal struc- 
ture, could they be derived on the assumption that the 
object is a point particle? The answer is in the affir- 
mative, provided that one is willing to introduce addi- 
tional assumptions and tolerate fields that are singular 
on the world line. The derivation sketched below was 
first presented by Mino, Tanaka, and Sasaki Q and then 
by Quinn and Wald 4]; the approach described here re- 
lies heavily on concepts and techniques introduced by 
Detweiler and Whiting [ToT ]. 



As we have seen in Sec. IIIII the gravitational pertur- 
bation produced by a point particle of mass m is ob- 
tained from the potentials of Eq. (JSJ by trace reversal: 
h a (3 — ip a ~ 5Sa/3(s 7l5 Vv)- I n Sec. Eq. © was as- 
sumed to hold in the external zone only [r > r e 3> m), 
but we now accept its global validity. We also accept the 
fact that the perturbation is singular on the world line 
(r = 0) . And we choose not to be bothered with the fact 
that while Eq. (S) was obtained by linearizing the Ein- 
stein field equations about the background metric g a p, 
the perturbation is obviously not small when r is smaller 
than m. 

Our first additional assumption is that the particle 
will move on a geodesic of the metric g a p — g a p + h a p. 
Writing down the geodesic equation in terms of tensorial 
quantities that refer to the background spacetime (such 
as the velocity vector u 11 = dz^/dr, which is normalized 
in the metric gap) returns the equations of motion 

= -^{gT + ^u v ) {2h„ X ; P ~ h Xp , v )u x uP. (17) 

These are superficially similar to Eq. I|15[) . but the differ- 
ence is important. While Eq. (|15fl involves tensor fields 
that are smooth on the world line, Eq. (|17fl involves 
highly singular quantities. These equations of motion 
are therefore meaningless as they stand, and h a p must 
be regularized before Eq. (|17fl can be evaluated. 

The regularization of a retarded field near its pointlikc 
source is a problem that has a long history. This issue 
was encountered by Dirac 0] m the context of the elec- 
trodynamics of a point electric charge in flat spacetime. 
By appealing to energy-momentum conservation across a 
world tube surrounding the charge's world line, Dirac dis- 
covered that regularization could be achieved by decom- 
posing the electromagnetic field into singular-symmetric 
"S" and regular-radiative "R" pieces; the "S" field would 
simply be removed from the retarded field and only the 
remainder — the "R" field — would be allowed to act on 
the particle. Dirac further discovered that in flat space- 
time, the "S" field is given by the following combination 
of retarded and advanced solutions to Maxwell's equa- 
tions: F%p = \{F^ + F*f). The "R" field, on the other- 
hand, is given by F^ = \{F^-F^). Because the "S" 
field satisfies the same field equations as the retarded field 
(with a singular source term on the right-hand side), it 
is just as singular as the retarded field on the world line; 
removing it from the retarded field produces a smooth 
field. This is confirmed by the fact that the "R" field 
satisfies the sourcefree Maxwell equations. 

Dirac's analysis was generalized to curved spacetime 
by DeWitt and Brehme but the proper decompo- 
sition of the retarded electromagnetic field into "S" and 
"R" parts was identified only recently by Detweiler and 
Whiting 0] . For reasons of causality, the flat-spacetime 
prescription (a superposition of half retarded and half 
advanced fields) does not work in curved spacetime: The 
retarded field depends on the particle's entire past his- 
tory, the advanced field depends on its future history, and 
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FIG. 1: Causal properties of the retarded and advanced solu- 
tions. The retarded field at a point x depends on the particle's 
history prior to (and including) the retarded point z(u). The 
advanced solution depends on the particle's history after (and 
including) the advanced point z(v). 

a linear superposition would depend on the full history, 
both past and future (see Fig. 1). This would give rise to 
equations of motion with unacceptable causal properties. 

The correct curved-spacetime prescription (identified 
by Detweiler and Whiting, and applied to the gravita- 
tional case) is to remove the "S" potential 

ips fi (x) = ^U af ! i , sl {x,x')ui u 5 

+ ^U a ^ s „(x,x")^"u s " 

-2m V'%{x,»)vf l u v dr (18) 

J u 

from the retarded potential. Here x' = z(u) stands for 
the retarded point on the world line associated with the 
field point x (and tensors with primed indices are eval- 
uated at that point), and x" = z(v) stands for the ad- 
vanced point (with a similar meaning for doubly-primed 
indices); within the integral z = z(t) stands for an ar- 
bitrary point on the world line, and the integral extends 
from the retarded point to the advanced point. Detweiler 
and Whiting have shown that 4>s^ ( x ) satisfies Eq. © - 
the same wave equation as the retarded potential — and 
comparison between Eqs. © and (|T%1) reveals that the 
"S" potential is just as singular as the retarded potential 
on the world line. The "R" potential is then defined by 

^(x) = r (x)-^f(x), (19) 

and it satisfies the homogeneous form of Eq. © ; this is 
smooth tensor field on the world line. The causal prop- 
erties of the "S" and "R" fields are illustrated in Fig. 2. 

Having introduced this unique decomposition of the 
retarded field into "S" and "R" pieces, Detweiler and 



FIG. 2: Causal properties of the "S" and "R" fields. The "S" 
field at a point x depends on the particle's history between 
(and including) the retarded point z(u) and the advanced 
point z(v). The "R" field depends on the particle's history 
up to the advanced point z(v). When x is brought to the 
world line the "R" field stays smooth and z(v) coincides with 
x; the field then depends on the particle's past history only. 

Whiting postulate that the "S" field exerts no force on 
the particle; the entire self-force arises from the action 
of the "R " field. This axiom can be seen as another un- 
avoidable assumption that must be introduced in order 
to make sense of the motion of point particles. Alterna- 
tively, it can be motivated by showing that the average 
of the "S" field on a spherical surface (as seen in the par- 
ticle's instantaneous rest frame) surrounding the particle 
is zero. (This calculation is performed in LRR yj.) In 
any event, the postulated equations of motion are 

^ = -^GT + uW) (2*&;p - h\,^u", (20) 

and one finds that the "R" field reduces to 

C;A = -4ro(«(A)M5 + u '" € + (21) 

on the world line, where h*"^\ is the tail field of Eq. i|16[l . 
Substituting Eq. (|21|l into Eq. I|20|l eliminates the terms 
involving the Riemann tensor, and we end up with the 
MiSaTaQuWa equations of Eq. JTSJ . We conclude that 
the equations of motion can indeed be derived on the 
basis of a point mass, at the price of introducing the 
geodesic postulate (which was not needed for a black 
hole) and the Detweiler- Whiting axiom (which also was 
not needed). 

Equation (|2U|) comes with a compelling interpretation. 
It states that the point particle moves on a geodesic of 
the metric g a p + h^p, which is smooth on the world line. 
Furthermore, because the "R" potential satisfies the ho- 
mogeneous form of Eq. ©, this metric is everywhere a 
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vacuum solution to the Einstein field equations. (Recall 
that g a p is itself a vacuum solution.) We therefore have 
a point mass moving on a geodesic of a well-defined vac- 
uum spacetime. 



V. NEWTONIAN SELF-FORCE 

The notion of a self-force exists also in the simple set- 
ting of Newtonian theory. And the Newtonian poten- 
tial also can be decomposed into "S" and "R" pieces. 
This pedagogical illustration was first presented to me 
by Steve Detweiler. The following presentation draws 
heavily from Ref. [l3| . 

Consider, in Newtonian theory, a large mass M at po- 
sition p(t) relative to the centre of mass, and a small 
mass m at position R(t). We assume that m -C M and 
the centre of mass condition reads mR + Mp = 0. We 
denote the position of an arbitrary field point by x, and 
r = \x\ is its distance from the centre of mass. We shall 
also let R = \R\ and p = \p\. 

We begin with a test-mass description of the situation, 
according to which the smaller mass moves in the grav- 
itational field of the larger mass, which is placed at the 
origin of the coordinate system. The background New- 
tonian potential is 



, . M 

<P (x) = 

r 



(22) 



and the background gravitational field is go = — V$ — 
—Mx/r 3 . In this description, the smaller mass m moves 
according to d 2 R/dt 2 = go(x = R). If the motion is cir- 
cular, then m possesses a uniform angular velocity given 
by flo 2 = M/R 3 , where R is the orbital radius. These re- 
sults are in close analogy with a relativistic description in 
which the smaller mass is taken to move on a geodesic of 
the background spacetime, in a test-mass approximation. 

We next improve our description by incorporating the 
gravitational effects produced by the smaller mass. The 
exact Newtonian potential is 



71/ 



\x - p\ 



\x-R\' 



(23) 



and for m <C M this can be expressed as 
(5$ (x), with a perturbation given by 



*o(a:)- 



6$(x) = 



M 



x-p\ 



M 
r 



\x-R\ 



(24) 



This gives rise to a field perturbation Sg that exerts a 
force on the smaller mass. This is the particle's "bare" 
self-acceleration, and the correspondence with the rela- 
tivistic problem is clear. 

An examination of Eq. (|24|l reveals that the last term 
on the right-hand side diverges at the position of the 
smaller mass. But since the gravitational field produced 
by this term is isotropic around R(t), we know that this 
field will exert no force on the particle. We conclude that 



the last term can be identified with the singular "S" part 
of the perturbation, 



$ s (x) 



\x-R\ 



(25) 



and that the remainder makes up the regular "R" poten- 
tial, 



$ R (£C) = -- 



M M 



\x — p\ r 



(26) 



The full perturbation is then given by <5<I>(a;) = $s( a; ) + 
$r(:e), and only the "R" potential affects the motion of 
the smaller mass. Once more the correspondence with 
the relativistic problem is clear. 

It is easy to check that to first order in m/M, Eq. I|2t)|) 
simplifies to 



R x 



(27) 



this simplification occurs because of the centre-of-mass 
condition, which implies that p is formally of order 
m/M <C 1. The "R" part of the field perturbation is 
then 



9r{x) 



m- 



3(R-x)x-r 2 R 



(28) 



and evaluating this at the particle's position yields a cor- 
rection to the background field go(x = R) = —MR/R 3 
given by gn(x = R) = 2mR/R 3 ; the force still points in 
the radial direction but the active mass has been shifted 
from M to M — 2m. For circular motion the angular ve- 
locity becomes O 2 = (M — 2m)/R 3 . This can be cast in a 
more recognizable form if we express the angular velocity 
in terms of the total separation s = R + p = (1 + m/M)R 
between the two masses. To first order in m/M we ob- 
tain f2 2 = [M + m)/s 3 , which is just the usual form of 
Kepler's third law. The "R" part of the field perturba- 
tion is therefore responsible for the finite-mass correction 
to the angular velocity. 

Notice that the "R" potential of Eq. (|27|l has a pure 
dipolar form. By contrast, the potentials (5$ and $s 
contain an infinite number of multipole moments. This 
observation should be kept in mind as the reader proceeds 
through the remaining sections of this contribution. 



VI. CONCRETE EVALUATION OF THE 
SELF-FORCE 

I turn next to the practical issues involved in a concrete 
evaluation of the gravitational self- force — the right-hand 
side of Eq. (|20|) . which shall be denoted a M [/i R ] — for a 
particle moving in the field of a Schwarzschild or Kerr 
black hole. 

The first sequence of steps are concerned with the com- 
putation of the (retarded) metric perturbation h a p pro- 
duced by a point particle moving on a specified geodesic 
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of the Kerr spacetime. A method for doing this was elab- 
orated by Lousto and Whiting 0] and Ori [l^ , building 
on the pioneering work of Teukolsky ^| , Chrzanowski 
|l7| , and Wald • The procedure consists of (i) solving 
the Teukolsky equation for one of the Newman-Penrose 
quantities tpo and '04 (which are complex components of 
the Weyl tensor) produced by the point particle; (ii) ob- 
taining from ipo or 04 a related (Hertz) potential '5 by 
integrating an ordinary differential equation; (iii) apply- 
ing to \& a number of differential operators to obtain the 
metric perturbation in a radiation gauge that differs from 
the Lorenz gauge of Eq. (J3J) ; and (iv) performing a gauge 
transformation from the radiation gauge to the Lorenz 
gauge. For a Schwarzschild black hole one can rely in- 
stead on the formalism of metric perturbations 0, |2j| , 
and the procedure simplifies. 

It is well known that the Teukolsky equation separates 
when 0o or 04 is expressed as a multipole expansion, 
summing over modes with (spheroidal-harmonic) indices 
I and m. In fact, the procedure outlined above relies 
heavily on this mode decomposition, and the metric per- 
turbation returned at the end of the procedure is also 
expressed as a sum over modes: 



K/3 = ^2 h U 



(29) 



(For each I, m ranges from —I to I, and summation of 
m over this range is henceforth understood. Indices on 
the metric perturbation and other tensors will now be 
omitted for ease of notation.) From the modes /i£, mode 
contributions to the self-acceleration can be computed: 
a[hg\ is obtained from Eq. l(2l7|) by substituting he in place 
of h R . These mode contributions are finite on the world 
line, but a[hi] is discontinuous at the radial position of 
the particle. The sum over modes does not converge, 
because the "bare" acceleration (constructed from the 
retarded field h) is formally infinite. 

The next sequence of steps is concerned with the reg- 
ularization of each a[h{\ by removing the contribution 
from h s . Starting with Eq. (|I8|I . the singular field can be 
constructed locally in a neighbourhood of the particle, 
and then decomposed into modes of multipole order £. 
This gives rise to modes a [/if] for the singular part of the 
self-acceleration; these are also finite and discontinuous, 
and their sum over £ also diverges. But the true modes 
a[hf] = a[hg] — a[hf] of the self-acceleration are continu- 
ous at the radial position of the particle, and their sum 
does converge to the particle's acceleration. 

The general structure of a [/if], as worked out in 
Refs. H IHHl HH13, is given by 

a" [/if] = (£+ + B» + 



convergent terms. 

(30) 

The "regularization parameters" A 11 , B 11 , and C 7 ' are in- 
dependent of I but they depend on the details of the 
spacetime and on the particle's state of motion. It is 
evident that the term involving in Eq. I|3U|I will di- 
verge quadratically when summed over £, that the term 



involving will diverge linearly, and that the term in- 
volving will diverge logarithmically; all other terms 
converge, and keeping additional terms produces faster 
convergence of the sum Y]( a [h>g]- 

The structure of Eq. H3U|) can easily be recovered in the 
Newtonian setting of the preceding section. The multi- 
polar decomposition of the singular potential of Eq. I|25|) 
is given by 



(31) 



where r < = min(r, R), 



max(r, R), n = x/r, R 



R/R, and Pi(n) are Legendre polynomials. Taking the 
gradient of Eq. I|31|) and then the limit x — > R returns 



9s 



m 
R 2 



E 



(32) 



for the radial component of the singular field (the an- 
gular components all vanish). This has the same form 
as Eq. (J3DJ, with A r = ^fm/R 2 , B r = -m/(2R 2 ), and 
C r = 0. The choice of sign in front of the {£ + |) term 
in Eq. Il.'i2[l comes from the two ways in which the limit 
r — > R can be taken: The upper (minus) sign corre- 
sponds to r — > R + while the lower (plus) sign corre- 
sponds to r — > R~ . This Newtonian calculation re- 
produces the leading-order term in a post-Newtonian 
expansion of the relativistic regularization parameters 

[Hmmmmmiii. 

The self-acceleration is obtained by first computing 
a[h(\ from the retarded metric perturbation, then com- 
puting the counterterms a[hf] by mode-decomposing the 
singular field, and finally summing over all a[hf] = 
a[h(] — a[hf\. This procedure is lengthy and involved, and 
thus far it has not been brought to completion, except 
for the special case of a particle falling radially toward 
a nonrotating black hole |28| . Another evaluation of the 
self-force, which did not rely on a mode decomposition, 
was carried out in the context of weak fields and slow mo- 
tions [2!|; this reproduced the standard post-Newtonian 
result. 

The procedure described above is lengthy and involved, 
but it is also incomplete when the background spacetime 
is that of a Kerr black hole. The reason is that the metric 
perturbations h l a g that can be recovered from 0o or 04 do 
not by themselves sum up to the complete gravitational 
perturbation produced by the moving particle. Missing 
are the perturbations derived from the other Newman- 
Penrose quantities: 0i, ip2, and ^3. While ipi and 1P3 
can always be set to zero by an appropriate choice of null 
tetrad, 02 contains such important physical information 
as the shifts in mass and ang ular- momentum parameters 
produced by the particle y(j . Because the mode decom- 
positions of ipo and 04 start at / = 2, one might say that 
the "£ = and £ = 1" modes of the metric perturbations 
are missing. It is not currently known how the proce- 
dure can be completed so as to incorporate all modes of 
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the metric perturbations. Specializing to a Schwarzschild 
spacetime eliminates this difficulty, and in this context 
the low multipole modes have been studied for the spe- 
cial case of circular orbits 0, . In view of the fact 
that the Newtonian self-force is purely 1 = 1 (refer back 
to the last paragraph of Sec.EJ, it is clear that these low 
multipoles cannot be ignored. 

Taking into account these many difficulties (and I 
choose to stay silent on others, for example, the issue 
of relating metric perturbations in different gauges when 
the gauge transformation is singular on the world line), 
it is perhaps not too surprising that such a small num- 
ber of concrete calculations have been presented to date. 
But progress in dealing with these difficulties has been 
steady, and the situation should change dramatically in 
the next few years. 



VII. BEYOND THE SELF-FORCE 

The successful computation of the gravitational self- 
force is not the end of the story. After the difficulties 
reviewed in the preceding section have all been dealt 
with and the motion of the small body is finally cal- 
culated to order m, it will still be necessary to obtain 
gauge-invariant information associated with the body's 
corrected motion. Because the MiSaTaQuWa equations 
are not by themselves gauge-invariant [they are formu- 
lated in the Lorenz gauge of Eq. ©], this step will ne- 
cessitate going beyond the self-force. 

To see how this might be done, imagine that the small 
body is a pulsar, and that it emits light pulses at regular 
proper-time intervals. The motion of the pulsar around 
the central black hole modulates the pulse frequencies as 
measured at infinity, and information about the body's 
corrected motion is encoded in the times-of-arrival of the 
pulses. Because these can be measured directly by a dis- 
tant observer, they clearly constitute gauge-invariant in- 
formation. But the times-of-arrival are determined not 
only by the pulsar's motion, but also by the propagation 
of radiation in the perturbed spacetime. This example 
shows that to obtain gauge-invariant information, one 
must properly combine the MiSaTaQuWa equations of 
motion with the metric perturbation. 

In the astrophysical context of the Laser Interferometer 
Space Antenna, reviewed in Sec. [HJ the relevant observ- 
able is the instrument's response to a gravitational wave, 
which is determined by gauge-invariant waveforms, /i+ 
and h x . To calculate these is the ultimate goal of the 
Capra program, and the challenges that lie ahead go well 
beyond what I have described thus far. To obtain the 
waveforms it will be necessary to solve the Einstein field 
equations to second order in perturbation theory. 

To understand this, consider first the formulation of 
the first-order problem. Schematically, one introduces a 
perturbation h that satisfies a wave equation Oh = T[z] 
in the background spacetime. Here T[z] is the stress- 
energy tensor of the moving body, which is a functional 



of the world line z(t). In first-order perturbation theory, 
the stress-energy tensor must be conserved in the back- 
ground spacetime, and z(t) must describe a geodesic. It 
follows that in first-order perturbation theory, the wave- 
forms constructed from the perturbation h contain no 
information about the body's corrected motion. 

The first-order perturbation, however, can be used to 
correct the motion, which is now described by the world 
line z(t) + 8z(t). In a naive implementation of the self- 
force, one would now re-solve the wave equation with a 
corrected stress-energy tensor, Oh = T[z + Sz], and the 
new waveforms constructed from h would then incorpo- 
rate information about the corrected motion. This imple- 
mentation is naive because this information would not be 
gauge- invariant. In fact, to be consistent one would have 
to include all second-order terms in the wave equation, 
not just the ones that come from the corrected motion. 
Schematically, the new wave equation would have the 
form of Oh = (1 + h)T[z + 5z] + (V/i) 2 , and this is much 
more difficult to solve than the naive problem (if only be- 
cause the source term is now much more singular than the 
distributional singularity contained in the stress-energy 
tensor). But provided one can find a way to make this 
second-order problem well posed, and provided one can 
solve it (or at least the relevant part of it), the wave- 
forms constructed from the second-order perturbation h 
will be gauge invariant. In this way, information about 
the body's corrected motion will have properly been in- 
corporated into the gravitational waveforms. 



VIII. CONCLUSION 

There has been significant progress toward solving the 
self-force problem over the last several years. The foun- 
dations (reviewed in Sees. II I II and IIV|> are now solid, 
thanks to the work of Mino, Sasaki, and Tanaka 01, 
Quinn and Wald and Detweiler and Whiting jlO| . 
The regularization parameters B^, and C M (intro- 
duced in Sec. IVI|I have been calculated for arbitrary mo- 
tion in the Schwarzschild and Kerr spacetimes by Barack 
and Ori 123, [2^, HSl and Mino, Nakano, and Sasaki 
|2lL l25l . and their expressions were independently veri- 
fied |2(1 127| . And finally, the reconstruction of the metric 
perturbation from the Teukolsky function (described in 
Sec. lVI|) is now better understood, thanks to the work of 
Lousto and Whiting |14j and Ori |15j . 

While this progress is encouraging, the challenges that 
lie ahead are still numerous. Among the outstanding is- 
sues are (i) the question of defining and calculating the 
u £ = and I = 1" contributions to the self-force in the 
Kerr spacetime (as was discussed in Sec. IVIj) . (ii) the 
design of useful gauge-invariant quantities that could be 
computed from the self-force and the metric perturba- 
tion (as was discussed in Sec. lVIIf) . and (iii) the incorpo- 
ration of the self-force into a consistent wave-generation 
formalism (as was also discussed in Sec. IVIl|) . These is- 
sues are fascinating, and in the next few years they will 
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be vigourously pursued by the Capra posse. 
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